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Abstract 

' We study genuinely exotic mesons with double charm and bottom flavor. We take a hadronic 

f-^ ; picture by considering the molecular states composed of a pair of heavy mesons, such as DD, DD* 

r"! , and D*D* for charm flavor, and BB, BB* and B*B* for bottom flavor. The interactions between 

heavy mesons are derived from the heavy quark effective theory. All molecular states are classified 
by I{J^) quantum numbers, and are systematically studied up to the total angular momentum 
J < 2. For the two-body problem, we numerically solve the coupled channel Schrodinger equations. 
We find the bound and/or resonant states for many quantum numbers. Some of them are new 
states which are not predicted by the tetraquark picture. 

PACS numbers: 12.39.Jh, 13.30.Eg, 14.20.-c, 12.39.Hg 
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I. INTRODUCTION 



Exotic hadrons which include multiquark configurations provide us with important infor- 
mation in hadron physics. There is a key to understand one of the most important problems 
in hadron physics; what are the constituent particles of hadrons, and what are the interac- 
tions among the constituent particles at relevant low energies. Nowadays the exotic hadrons 
are studied, not only in light fiavor sectors, but also in heavy flavor sectors with charm and 
bottom quarks l|-l5|. Recent experimental observations of heavy exotic hadrons, Dgj, X, Y, 
in charm sector, and Yb, in bottom sector have motivated many physicists toward 
the question about the possible new dynamics in the heavy hadrons. Many of those hadrons 
have masses, decay widths and branching ratios, which may not be explained as normal 
hadrons, such as qq and qqq. Especially Z^ and Z^ which have an isospin one cannot be 
assigned as simple quarkonia, but there should be additional up and down quark degrees 
of freedom. Their quark contents will be considered as QQqq with a heavy quark Q and a 
light quark q, because qq can form isospin one. 

In the present paper, we consider a new hadron state Tqq whose quark content is QQqq 
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39| . Tqq is a system containing two heavy quarks and is genuinely exotic, be- 



cause there is no annihilation process of a quark and antiquark pair. As discussed in many 
literatures j6-29|, Tqq is considered as a tetraquark state, in which the effective degrees of 
freedom are constituent quarks. There, it is shown that Tqq can be a stable object due to 
the strong attraction in qq which may form a stable scalar diquark 30l432 L Tqq may be a 



deeply bound state, which does not decay through strong interaction 22-2^. The study of 
Tqq is also useful to understand the color superconductivity in high density quark matter, 



because the diquark correlations are important to cause a diquark condensate 
far tetraquark states such as T^^, T]^ and T^j^ states with I(J^) = 0(1 



33|-|35|. So 
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as well as T^^ with 



24|. 



I{J^) = 0(0+) have been discussed as stable objects. 

However, the deeply bound tetraquark is not the only possible picture for Tqq. When four 
quarks are present, they may form Qq hadronic clusters which alternatively become relevant 
degrees of freedom instead of diquarks. This is the basis of hadronic molecular states which 
are expected to appear near the hadronic threshold region 37H39j. In the present Tqq 



system, two mesons composed by Qq, the pseudoscalar meson P ~ (Qq)spino and the vector 
meson P* ~ (Qq)spini, can become effective degrees of freedom as constituents. In the heavy 
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quark limit, the pseudoscalar meson P and the vector meson P* become degenerate in mass, 
and hence both of them should be considered. Hereafter we introduce the notation P*^*^ for 
P and P*. 

The interaction between two P^*^'s is supplied by the meson exchange interaction. Because 
P^*) carry a finite isospin, the pion (vr) exchange interaction works dominantly at long 
distances. The existence of a pion is a robust consequence of spontaneous breaking of chiral 



symmetry [36|. Thus, the mass degeneracy of P and P* is concerned with the heavy quark 

(spin) symmetry, and the interaction between two P'^*^'s with chiral symmetry. Therefore, 

those two symmetries are important when Tqq is considered as the hadronic molecule of 
p(*)p(*)_ 

Until now there have been several studies for hadronic bound states of p(*)p(*). The 
pioneering work by Manohar and Wise showed that the PP(*) mesons could form a bound 
state [37|, and the result has been followed by other authors |38|, l39|. However, the quantum 



numbers I{J^) which have been discussed so far are limited. Moreover, resonant states of 
p(*)p(*) have not been studied so far. In the present work, we investigate systematically 
various states with possible quantum numbers I{J^) with J < 2, and discuss both bound 
and resonant states of p(*)p(*) by considering the fully coupled channel problem. 

This paper is organized as follows. In section 2, we give the interaction between two P'-*^ 
mesons with respecting the heavy quark symmetry and chiral symmetry. We introduce the 
two types of potentials, the vr exchange potential and irpuo exchange potential. In section 
3, we apply these potentials to the p(*)p(*) systems and solve the Schrodinger equations 
numerically to search the bound and/or resonant states. In section 4, we compare our 
results from the hadronic molecule picture with the previous results from the tetraquark 
picture. We summarize our discussions in the final section. 



II. INTERACTION WITH TWO MESONS WITH DOUBLY HEAVY FLAVOR 

The dynamics of the hadronic molecule of p(*)p(*) respects two important symmetries; 
the heavy quark symmetry and chiral symmetry. The heavy quark symmetry induces the 
mass degeneracy of P and P* in the heavy quark limit. Because of this, we have to consider 
the channels of degenerate pairs, such as PP, PP*, P*P and p*p*^ leading to the mixing of 
them; PP*-P*P, p*p*-P*p*, PP-P*P*, pp*.p*p*. 
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As for the interaction, the one pion exchange potential (OPEP) exists between P(*)'s at 
long distances. The OPEP is provided by the PP*7r and P*P*7r vertices whose coupling 
strengths are equally weighted thanks to the heavy quark symmetry. We note that there is 
no PPtt vertex because of the parity conservation. Therefore, at long distances, the OPEP 
plays an important role to realize the channel mixings in p(*)p(*). At short distances, heavier 
mesons also should be considered. 

To derive the p(*)p(*) potential, we em ploy the effective Lagrangians based on the heavy 



quark symmetry and chiral symmetry 
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45[ . They describe the interaction between heavy 



mesons which are given by the exchange of pion and vector mesons {v = p,uj). The interac- 
tion Lagrangians are given as 

^.pp- = 2f (Pin; + p:iPb)d^n^b , (1) 

Jit 

C^P^p, = 2i^e^^^^v^P:lP:^d,Kb: (2) 

Jvr 

C,vPP = -V2f3gvPbPlv ■ pba, (3) 
C^PP, = -2V2\gyv,e^''-^ {pIp;^ - P^P,) d^ba , (4) 

+l2^\gyP:lP;^{d^{nba - d'^mba) , (5) 

where P = D, B and P* = £)*, 5*. The subscripts a and b are for light flavor indices, up 
and down, and is a four- velocity which will be fixed as = (1, 0) below. The pion and 
vector meson (p and u) fields are defined by 




Following Ref. 45|], the coupling constants in the interaction Lagrangians are given as 



g = 0.59, P = 0.9, A = 0.56 GeV-\ U = 132 MeV, gv = (7) 
where f.^ is the pion decay constant and rUp is the mass of the p meson. The OPEPs are 
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derived by the interaction Lagrangians ([T]) and ([2]) as follows: 



Vp*P*^P*P* 



V, 



PiP2^PiP^ 



Jtt 



[^i -£2 "^7r) + 5'e*,£. T(r; m^)] frfa, 



1 



(8) 
(9) 
(10) 

(11) 



where m,r is the pion mass. Here three polarizations are introduced for P* as defined by 
e^^) = (=f1/-\/2, ±z/v^, 0) and = (0,0,1), and the spin-one operator T is defined by 
Tl^ix = is'^-^'^s'^^^'^ e'^^\ As a convention, we assign e^^^ for an incoming vector particle and 
^W* for an outgoing vector particle. Here fi and 72 are isospin operators for and P2*''; 
^1 ■ ^2 = ~3 and 1 for the J = and 1=1 channels, respectively. We define the tensor 
operators 

5',*,,, = 3(£(^^)*-f)(£(^2)-f) 
Sti,T2 = 3(fi-f)(f2-f) - fi-fs. 



(12) 
(13) 
(14) 



(15) 



where f = r/r is a unit vector between the two mesons. 

The p meson exchange potentials are similarly obtained from the interaction Lagrangians 

(ED-dSD, 



^PlP2^PlP2 



^PlP,*^PlP,* 



(3gv 

13 9v 
2m„ 



-C(r;mp)fi-f2, 



-C{r;mp)fi-f2, 



Vp^p;^p*P2 = (2A^v)% [2e*-e2C{r-mp)-Sei,s2T{r;mp)] n-rs. 



^1 ^2 "'■^l ^2 



'^PlP2^Pl*P2* 
yp^P*^P*P* 



(2A5fy) - 2TrT2C{r;mp)-ST,,T2T{r;mp) 



ri-T2 



+ 



^9v_ 
2mn 



-C{r;mp)fi-f2, 



{2\gvf-[2ele;C{r-mp)-Se,^,e,J{r-mp)\ fi-ra. 



{2\g 



V 



2e*-T2 C{r; mp)-Sei,T2 T{r; nip) 
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Ti-r2. 



(16) 

(17) 
(18) 

(19) 

(20) 
(21) 



The oj meson exchange potentials are obtained by replacing the mass of p meson with the 
one of a; meson and by removing the isospin factor Ti-f2. The OPEP's of p(*)p(*) differ from 
the ones of p(*)p(*) in that the overall signs are changed due to G-parity. The situation is 
the same with uj meson exchange potentials, while p meson exchange potentials of p(*)p(*) 



are not changed because the G-parity is even 
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In the above equations, C{r;mh) and T{r\mh) are defined as 

C{r-mn)= I ^^^L^e^^-F(g;m,), (22) 
J {2'kY + ml 

T{r- m^)S,,{r) = [ ^ S^Me'^'^^m mn), (23) 

with 5*12 (x) =3{(Ti-x){a2-x) — (j\-a2- We introduce the monopole type form factor at each 
vertex to take into account of the size effect of P*^*-* mesons. So the function reflected form 
factors is defined as 

m^.)=(^%f^S . (24) 



where m/j and q are the mass and three-momentum of the exchanged meson h (= 'K,p,u) 
and Ap is the cut-off parameter. The cut-off parameter Ap are determined from the size 



of P estimated from the constituent quark model as discussed in Refs. 46 



. The cut-off 



parameters are Ap, = 1121 MeV and Ap = 1070 MeV when the vr exchange potential is 
employed, while Ap = 1142 MeV and Ap = 1091 MeV when the vrpa; is employed. 



III. BOUND AND RESONANT STATES 

Let us classify all the possible quantum numbers of the p(*)p(*) systems with isospin 
/, total angular momentum J (J < 2), and parity P. We show the quantum numbers 
I{J^) and the channels in the wave functions in Table [B It is noted that the wave func- 
tions must be symmetric under the exchange of the two P*^*^ mesons. We use the notation 
"^^^^L J to indicate the states with the internal spins S and angular momentum L. For ex- 
ample, the I{J^) = 0(1+) state is a superposition of four channels; ^ (PP* — P*P) (^^i), 

(PP* - P*P) i^Di), F*F*{^Si) and F*F*i^Di). It is important to consider the mixing of 
all possible channels. Especially the tensor force plays a crucial role to mix the channels 
with different angular momenta, L and L ± 2. We obtain the Hamiltonian in a matrix 
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form of these coupled channels. The explicit forms of the Hamiltonian are summarized in 
Appendix A. 

We numerically solve the coupled channel Schrodinger equations for the states with each 



quantum number. As a numerical method, the renormalized Numerov method [50| is adopted 
to solve the coupled second-order differential equations. The resonant states are identified 
by the behavior of the phase shift 5 as a function of the scattering energy E. The resonance 
position Er is defined by an inflection point of the phase shift 6{E) and the resonance width 



by Tr = 2 / {dS / dE) E=Er following Ref. j47|, |48|, |5l| . To check the consistency of our numerical 



calculations, we also adopt the complex scaling method (CSM), in which the resonant state 



is defined as a pole in the complex energy plane (5l|, |52| . We obtain an agreement in the 
results of the renormalized Nemerov method and the CSM. 

We summarize our numerical results for D^*^D'^*) bound/resonant states in Table HTl and 
FiglH In D*^*''D*^*) states, we find several bound and/or resonant states in J = 0, while there 
is no bound state in / = 1. In general, the attractive force of pion exchange in / = 1 is three 
times weaker than in / = due to the isospin factor. As a numerical result, D^*''D^*'' bound 
states in / = 1 are not obtained but only resonant states are. 

Let us look at our results one by one for each quantum number in detail. In the followings, 
most of the numerical values are shown for the case of the irpu potential, because the results 
from the npu potential are generally not so different from those from the vr potential, except 
for the 0(2^) state. 

0(0~) This state has only one channel and the pion exchange potential is attractive as shown 
in Eq. (IA14p . As a result, this state is a very deeply bound state with binding energy 
132.1 MeV measured from the DD* threshold. 

0(1"*") The pion exchange potential is repulsive for diagonal components as shown in 
Eq. flAlSp . However this state has four components and the mixing of the S- and 
Z^-waves causes the strong tensor attraction from the off-diagonal components of the 
potential. Consequently, there is a deeply bound state with binding energy 62.3 MeV 
measured from the DD* threshold. 

0(1") There are twin shape resonances with the resonance energy 17.8 MeV and the decay 
width 41.6 MeV for the first resonance, and the resonance energy 152.8 MeV and the 
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decay width 10.6 MeV for the second. The resonance energies are measured from the 
DD threshold. Those resonances are formed by the centrifugal barrier in the P-wave. 

0(2"'") This state contains only D-wave components. The potential is weakly attractive. 
Nevertheless, due to the centrifugal barrier in the D-wave, there is a resonant state at 
the energy 33.7 MeV from the DD* threshold, but the decay width 196.3 MeV is very 
wide. 

0(2"") When the OPEP is employed, there are twin resonant states with the resonance energy 
0.1 MeV and the decay width 0.02 MeV for the first resonance, and the resonance en- 
ergy 118.0 MeV and the decay width 23.4 MeV for the second from the DD* threshold. 
When the effects of p and u meson exchange are included, the first resonance becomes 
a weakly bound state with the binding energy 4.3 MeV, because the p meson exchange 
enhances the central force attraction of the pion exchange. The u meson exchange 
plays a minor role due to the isospin factor, although this contribution suppresses the 
attractive central force. The second resonant state with the resonance energy 112.1 
MeV and the decay width 26.6 MeV is not affected very much. 

1(0"") This is the only 1 = 1 state in D^*''D^*\ The interaction in J = 1 are either repulsive 
or only weakly attractive as already discussed. Nevertheless, due to the P-wave cen- 
trifugal barrier, we find twin shape resonances; the resonance energy 2.3 MeV and the 
decay width 37.4 MeV for the first resonance, and the resonance energy 144.2 MeV 
and the decay width 34.4 MeV for the second, from the DD* threshold. 

Next we discuss the results for the B*^*'B*^*'' states in Table IIIII and Figs. |2] and [31 In 
Figs. |2] and [31 only states around the threshold region are plotted. For B(*^B(*^ states, we 
use the same coupling constants, and change only the masses of heavy mesons with small 
difference of the cutoff parameters. At first glance B*^*''B*^*'' states have many bound and 
resonant states in comparison with d(*)d(*) states. There are two reasons. First, the kinetic 
term is suppressed in the Hamiltonian because the reduced mass becomes larger in bottom 
sector. Second, the effect of channel-couplings becomes more important in the bottom sector, 
because a pseudoscalar meson B and a vector meson B* become more degenerate thanks to 



the heavy quark symmetry; similar discussion has been done in Refs. 



46 



-48 



We note that, 



as a consequence of the strong attraction, several new states appear in the B(*)B(*) states 
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m isospm triplet; /(J^) = 1(0+), 1(1"), 1(2+) and 1(2"). The corresponding states are not 
obtained in the D(*)D(*) states. 

Here we comment on the similarities in the energy spectrum in the charm and bottom 
sectors. When we look at the / = states in Figs. [T] and [21 and the l(O^) states in Figs. [T] 
and [31 we find qualitatively the correspondence between the states in charm and bottom 
sectors. For instance, I{J^) = O(O^) has one state, and I{J^) = 0(1") has two states in 
both charm and bottom sectors. The correspondence between the two sectors may lead us 
to the idea that the energy spectrum around the thresholds does not depend on the heavy 
fiavor very much, when the PP, PP* and P*P* thresholds are close each other. If we look 
at their levels in the threshold region in more detail in / = 0, the states in the charm sector 
are those of principal quantum number n = 1 (lowest) or n = 2, while those in the bottom 
sector plotted in Figs [2] and [3] are with n = 3. The correspondence is more clearly seen in 
the 1(0") states in J = 1, where these states happen to be shallow bound states or resonant 
states both in charm and bottom sectors. Therefore, we observe that the states appearing 
around the thresholds exist as hadronic molecules and share similarity in their spectrum. 
The states with stronger binding far below the thresholds with smaller n may no longer 
survive as hadronic molecules, but rather tetraquarks, or mixtures of hadronic molecules 
and tetraquarks 5J]. The tetraquark will be briefiy discussed in the next section. 



IV. HADRONIC MOLECULES OR TETRAQUARKS? 

In the present paper, we have considered Tqq as a hadronic molecule composed by p(*)p(*) 
in which one pion exchange potential induces a dominant attraction. On the other hand, as 
mentioned in the introduction, Tqq can also be considered as a tetraquark QQqq in which a 
diquark qq provides a strong binding energy. We briefiy comment on the difference of those 
two pictures. 

The different features between the hadronic molecule and tetraquark pictures are seen 
in their masses. For hadronic molecules to be valid, hadron constituents are sufficiently far 
apart such that their identities as hadrons must be maintained. They can not be too close 
to overlap each other. Therefore, masses of hadronic molecules should appear around their 
threshold regions. In contrast, tetraquarks may be strongly bound and become compact 
objects as genuine quark objects. Thus, their natures are differentiated in the masses, small 
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TABLE I: Several channels of p(*)p{*)(25+i^^) each state with / and for J < 2. 



I 




components 




Q- 


1 Cpp* I p*p)('3p \ 




1+ 


1 /pp* _ p*p\ (3g \ 1 Cpp* _ p*p^ fSn N p*p*(3c N P*P*(3£) ) 





f- 


pp(lR) 1 (pp* +P*P) (3p \ P*P*(lPi), P*P*(5r) P*P*(5p ) 




2+ 


1 Cpp* _ p*p\ (3^ \ P*F*(^D2) 




2" 


1 /pp* ^ p*p\ (Sp y 1 /pp* _^ p*p\ /3p N p*p*(5p y p*p*(5p ) 




0+ 


PPf^S'n) P*P*f^S'n) P*P*f^i:>n) 




0~ 


1 /pp* _ p*pN /3p \ p*p*(3p ) 


1 


1+ 


1 /pp* _^ p*pN /3 o N 1 /pp* _^ p*pN (3^ ) P*P*(5Z)i) 




1" 


(PP* - P*P) P*P*(3P;^) 




2+ 


PP(iZ)2), ^ (PP* + P*P) (^£'2), P*P*(^L>2), P*P*(^52), P*P*(5L>2), P*P*(^G2) 




2" 


(PP* _ p*p) (3p2), ^ (PP* _ p*p) (3^2), P*P*(3P2), P*P*(3F2) 



binding energy of order ten MeV or less, or larger one. Although it is tempting to seek for a 
framework to cover both scales, such a problem is out of scope of the present paper. Instead, 
we compare the results from the two pictures and just clarify the differences between them. 

In the hadronic molecule picture by p(*)p(*)j as presented in the previous section, we 
obtain, not only the bound states, but also the resonant states in many I{J^) quantum 
numbers. Both in charm and bottom sectors, it is remarkable that there are even the twin 
resonances in several quantum numbers, 0(1"), 0(2") and 1(0"). In the tetraquark picture 
including a diquark model 23|, l2J|, in contrast, only two bound states rnliJ^) = 0(1"^) and 
1(2") have been predicted so far, as discussed for an example in Refs. 27H29| as recent works. 
For 0(1"*"), the predicted binding energy can be around 70 MeV from the DD* threshold, 
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TABLE II: The energies of D(*)d(*) states with I{J^) with J < 2. The energies E can be either 
pure real for bound states or complex for resonances. The real parts are measured from the 
thresholds as indicated in the second columns. The imaginary parts are half of the decay widths of 
the resonances, T/2. The values in the parentheses for the bound states are matter radii in units 
of fm. 



/ 






E [MeV] 


/I -pULcllLld,! 


•TT /~t / . 1 T~* /~»T on T 1 O 1 

/I , fJ^ U/-pULtJllLifii 







DD* 


oyj.^y ±.yj J 




1+ 


DD* 


—A^ 7(0 Q'] 


— fi9 '^(0 R') 


1" 


DD 


iq zL _ ift'^ 0/9/9 A — ?"^7 4/9 


1 7 S _ 741 /9 1 ^9 R — vl /9 


2+ 


DD* 


34.5 -«183. 1/2 


33.7- il96.3/2 


2^ 


DD* 


0.1 - iO.02/2, 118.0 - i23.4/2 


-4.3(1.6), 112.1 - i26.6/2 


1 


0+ 


DD 


no 


no 


0" 


DD* 


3.7-i27.7/2, 143.8 -^40.0/2 


2.3 - i37.4/2, 144.2 - i34.4/2 


1+ 


DD* 


no 


no 


1~ 


DD* 


no 


no 


2+ 


DD 


289.4 -ilO. 9/2 


no 


2- 


DD* 


no 


no 



which is accidentally similar to our values 62 MeV. For 1(2^), the predicted binding energy 
can be around 27 MeV. The 1(2") state cannot be found in our present study. Thus, 
we observe that the energy spectrum of the hadronic molecule picture by p(*)p(*) is quite 
different from that of the tetraquark picture. 
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TABLE III: The energies of B(*)b(*) states with I{J^) with J < 2. (Same convention as TablelUl) 



/ 




threshold 


E-iT/2 [MeV] 








TT-potential 


vr, yO, w-potential 




Q- 


BB* 


178.0(0.6), -32.0(1.2) 


-305.9(0.5), -77.5(0.9), -3.3(2.5) 




1+ 


BB* 


-179.7(0.5), -25.7(1.2) 


-201.5(0.5), -33.6(1.1) 





1" 


BB 


-125.5(0.6), -39.1(0.8), 
1.9 - i3.3/2, 52.8 - il2.8/2 


-164.4(0.5), -98.9(0.6), 
-3.1(1.6), 35.0 - ill.8/2 




2+ 


BB* 


-51.2(0.8), 5.5 -i5.9/2 


-60.6(0.8), 5.7- il3.2/2 




2- 


BB* 


-147.5(0.6), -68.7(0.8), 
7.6-i4.3/2, 26.9 -i20. 2/2 


-196.5(0.6), -84.1(0.7) 
0.5 - i8.5/2, 7.6 - i4.3/2, 26.9 - i20.2/2 




0+ 


BB 


-18.1(0.9), 


-33.9(0.5) 




0" 


BB* 


-50.5(0.8), 0.7 - i3.0/2, 46.7 - il.9/2 


-5.9(1.4), 38.5 - i4.6/2 


1 


1+ 


BB* 


-38.1(0.8) 


no 




1" 


BB* 


no 


11.7-«11.0/2 




2+ 


BB 


23.0 - i4.4/2 


62.4-^52.5/2 




2- 


BB* 


2.0 - i3.3/2, 63.7- i7.6/2 


2.3-i4.7/2 



V. SUMMARY 

We have discussed exotic mesons with double charm and bottom flavor whose quark 
content QQqq is genuinely exotic. We have taken the hadronic picture, and considered 
molecular states of two heavy mesons p(*'''s (a pseudoscalar meson P = D, B, and a vector 
meson P* = D*, B*). With respecting the heavy quark and chiral symmetries, we have 
constructed the it exchange potential and the irpu exchange potential between the two 
heavy mesons. To investigate the bound and/or resonant states, we have numerically solved 
the coupled channel Schrodinger equations for the p(*)p(*) states with I{J^) for J < 2. 

As results, we have found many bound and/or resonant states in both charm and bottom 
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D*D* 
(4017) 



4020 



3988 (34) 
(27) 



3909 

DD* ^ ^'''^ 3871 J87i 

(3876) (11) (37) 

3813 



3743 3743 

DD —T- 

(3734) (^^) 

1(F) 0(0) 0(1+) 0(1) 0(2+) 0(2) 1(0) 



FIG. 1: The D(*)d(*' bound and resonant states with I{J^). Solid lines are for our predictions 
for the energies E of the bound and resonant states, as denoted above the lines, and the values in 
parentheses below the lines denote the decay width T of the resonances when the vrpo; potential is 
employed. Mass values are shown in units of MeV. 

sectors. The D(*^D(*^ bound and resonant states have moderate energies and decay widths 
around the thresholds in several channels with quantum numbers; O(O^), 0(1"'"), O(l^), 
0(2+), 0(2~) and 1(0~). The B'-*^B(*' states have more bound and resonant states with 
various quantum numbers. Several new states appear in the b(*''B(*^ states in isotriplet 
states, such as 1(0"'"), 1(1~), 1(2"'") and 1(2~"), which cannot be found in the charm sector. 
By contrast to the D(*)D(*) states, some B^*^B'-*^ states are very compact objects with a large 
binding energy much below the thresholds. Perhaps, these states cannot survive as hadronic 
molecules and more consideration of quark dynamics such as tetraquarks is required. 

Experimental studies of those exotic hadrons should be performed in the coming fu- 



ture. The double charm production in accelerator facilities will help us to search them [53|. 
Recently it has been discussed that the quark-gluon plasma in the relativistic heavy ion col- 
lisions could produce much abundance of exotic hadrons including the exotic mesons with 



double charm l5a. 



56| . Those experimental studies will shed a light on the nature of the 
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B*B* -- - 

(10650) 

10631 
(20) 

10610 1Q612 
gg* 10601 (13). (4) 

(10604) 10594 

(12) 

10571 

BB 1Q556 
(10559) 

1(F) 0(0) 0(1+) 0(1) 0(2+) 0(2-) 



FIG. 2: The B(*)b(*) bound and resonant states around the thresholds with I{J^) in I = 0. 
(Same convention as Fig. [TJ) 



B*B* 
(10650) 10643 

(5) 

10621 

10616 — 

10607 

(lSV"""i^ (5)"" 

BB 

(10559) 

I(JP) 1(0) 1(1) 1(2+) 1(2) 



FIG. 3: The B'^*)B^*) bound and resonant states around the thresholds with I{J^) in I = 1. 
(Same convention as Fig. [TJ) 
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exotic mesons with double charm and bottom flavor, and provide hints to the fundamental 
questions of the strong interaction in hadron physics. 
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Appendix A: Hamiltonian 

The Hamiltonian is a sum of the kinetic term and potential term as. 



for the Tipuj potential. 

The kinetic terms with including the explicit breaking of the heavy quark symmetry by 



Hl[jP) = Kil^jP) + V[^jPy 



(Al) 



for the TT exchange potential only, and 




(A2) 
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the mass difference Ampp* = mp. — mp are 

i^'o(i+) = diag ( -7— ^ — Ao, -— ^ — A2, Aq + Ampp., -— ^ A2 + Ampp, 

y zmpp. impp. 2mp*p. zmp.p. 

(A4) 

^o(i-) = diag ( -— J— Ao, -— ^ — Ai + Ampp., 
y zmpp zmpp* 

1 1 1 \ 

Ai + 2Ampp*, r Ai + 2Ampp*, z A3 + 2Ampp. , (A5) 



2mp.p* 2mp.p. 2mp.p. 
Ko(2+) = diag ( -— ^A2, -— ^ A2 + Ampp. ) , (A6) 



1 


■A2,- 


1 


2mpp* 


2mp.i 


1 


■Ai,- 


1 


2mpp. 


2mpp 



Ko{2-) = diag ( -— ^ — Ai, -— ^ — A3, -7—^ Ai + Ampp., -— A3 + Ampp, 

\ 2mpp. 2mpp. 2mp.p. 2mp.p. 

(AT) 

iri(o+) = diag ( -— — Ao, -— Ao + 2Ampp. , -— A2 + 2Ampp. ) , (A8) 

y 2mpp 2mpp. 2mpp* y 

^■"*> = ^'^s {-^^"' -^^'^ -2^"^' + ^"-) • 

Ki{2+) = diag ( -7-^^ — A2, -—^ — A2 + Ampp., -7:-^^^ A2 + 2Ampp., 

y 2mpp 2mpp. zmp.p. 

-7—^ Ao + 2Ampp., -7— ^ A2 + 2Ampp., -7—^ A4 + 2Ampp. ) , (A12) 

2mp.p. 2mp.p. 2mp.p. J 

Ki(2-) = diag ( -77-3^^ — Ai, -7—^^ — A3, -77-^^^ Ai + Ampp., -7—^^ A3 + Ampp. 

V 2mpp. 2mpp. 2mp.p. 2mp.p. 



(A13) 



where A, = |^ + f ^ - with integer / > 0, 1/mpp = 1/mp + 1/mp, 1/mpp. 

1/mp + 1/mp., 1/mp.p. = 1/mp. + 1/mp.. 
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The TT exchange potentials for each I{J^) states are 



o(o-) 



"^0(1 + ) 



^0(1- 



'^0(2+) 



0(2- 



1(0+) 



y-rr 



^1(1+) 



yn 



^1(2+) 



yTT 

^1(2-) 



(Vb + 2Vt) , 

( -Vc V2Vt 2Vc V2Vt \ 

y/2VT -Vc-Vt v^Vt 2Vc - Vt 

2Vc V2Vt -Vc V2Vt 

V V2Vt 2Vc - Vt V2Vt -Vc - Vt J 

( 






-2^1 Vt 





Vc - Vt 


-3^1 Vt 
-sVIvt 



-VsVc 


-2Vc 



-3^|Vt 



^Vt Vc - |Vt 



iVT 



Vt 



S^fVT 

-3^1 Vt 

fVT 

^Vt 
Vc - f Vt ) 



-Vc + Vt 2Vc + Vt 
2Vc + Vt -Vc + Vt 



/ Vc + iVT 



3V& 



^Vt Vc + I 



Vt 
Vt 



^Vt 

5 



Vt 

— Vt ^Vt Vc + |Vt 



-^Vt \ 



'±Vt 



3V^ 



6 



Vt 



Vt -2^Vt 

-VsVc V6Vt 
VSVc -2Vc -V2Vt 
\^ V6Vt -V2Vt Vc - 2Vt ) 

-Vc- 2Vt 2Vc - 2Vt 
2Vc - 2Vt -Vc- 2Vt 

' Vc -V2Vt -V^Vt \ 
-y/2VT Vc + Vt -VsVt 
-VgVt -VsVt Vc-Vt ) 

-Vc + Vt -2Vc - Vt 
2Vc - Vt -Vc + Vt 



Vc 



Vt / 



(A14) 



(A15) 



(A16) 



(A17) 



(A18) 



(A19) 



(A20) 



(A21) 



(A22) 





Vc-Vt 







-\/3Vc 
Vt 

fVT 

ISVt 



( -Vc - \Vt ^Vt 



-V3Vc 


-2Vc 



S^fVT 
-yfVT 



IVt 



-^Vt 



35 



Vt \ 



%/7 



Vt 



,Vt 



3^1 Vt 


-^fVT 


Vc 













Vc + f Vt 


-^Vt 









^Vt 

7\/5 ^ 


Vc - f Vt / 



(A23) 



2Vc - iVT 



3a/6 



2Vc 
3v^ 



Vt 
|Vt 
Vt 



-Vc-I 

3v^ 



Vt 



2Vc - |Vt 



3v^ 
5 

-Vc - 



Vt 



^Vt \ 
2Vc - |Vt 
Vt 



3\/6 



^Vt 



(A24) 



-Vc - |Vt / 
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The p and 00 potentials are 



yv 



yv 

'^0(2+) 



yv 

"0(2") 



yv 

^(0+) 



yv 

"1(0-) 



yv 

S(l+) 



yv 
l(l-) 



yv 

S(2+) 



yv 

"1(2-) 



{2VS - + v^') 

( 



V 




4y(j + V5; 


/ 












2vg + + ^(5' 













3yf^T" 


\ 




-3^1^? 



-V2y^ 

-2Vg + Vg' 

-V2y^ 



-2v^y(^ 



4y-- 



1/ 5 '^T 





\ 
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\ 


3^1 Vt" 
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5 
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/ 
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5 
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5 

6V^ 
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' 5 
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5 
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5 "^T 
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\ -yey^' 
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y^ 
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-2y^ + 2y^ + y^' 
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where the central and tensor potentials are defined as, 

= (v^|-) \c{r-m^)n-r2, (A36) 

= (v^|-) ^T(r;m^)frf2, (A37) 

= {2\gyf^Cir;mp)n-f2, (A38) 

= {2Xgvf^Cir;m^), (A39) 

V4 = {2Xgvf^T{r;m,)n-T2, (A40) 

= {2\gyf^Tir;m^), (A41) 
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